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Abstract 

Let Aj be the J-th eigenvalue of Sturm-Liouville systems with separated boundary conditions, we 

build up the Hill-type formula, which represent n(l “ as a determinant of finite matrix. This is the 

J 

first attack on such a formula under non-periodic type boundary conditions. Consequently, we get the 

Krein-type trace formula based on the Hill-type formula, which express 2 as trace of finite matrices. 

j ‘ 

The trace formula can be used to estimate the conjugate point alone a geodesic in Riemannian manifold 
and to get some infinite sum identities. 
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1 Introduction 

In this paper, we will consider the eigenvalue problem for the Sturm-Liouville systems 

- {Py + Qy)' + Q^y + {R + m)y = 0, (i. i) 

where 2 is a continuous path of n x n matrices, and P,R,Ri are continuous paths of n x « symmetric matrices 
on [0, rj. Instead of Legendre convexity condition, we assume that for any t e [0, T], P{t) is invertible. 

The eigenvalue problem of the Sturm-Liouville systems depends on the boundary conditions. There 
are two important type boundary conditions, periodic type and separated type. Lor the literature in ?i-body 
problem, readers can refer to |l5l,|[8l. The eigenvalue problem for S-periodic boundary value problem, that 
is, y(0) = Sy{T) for some orthogonal matrix S , was studied in Q. In the present paper, we will consider the 
separated boundary conditions, which includes the homogenous Dirichlet, Neumann and Robin boundary 
conditions. More precisely, let Aq, Ai be two Lagrangian subspaces of (R^",a>o) which are phase spaces 
with standard symplectic structure. Set x = Py + Qy, z = (x,y)^, and the separated boundary condition is 
given by 


z(0)eAo, z{T)eAi. (1.2) 
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In order to understand the eigenvalue problem of the system (I1.1H1.2I) . for the eigenvalues Aj we will 
build up a formula with the form 0/(1 ~ = det(Al), where the matrix M depends mainly on the 

monodromy matrix and the Lagrangian subspaces Aq, Ai. We called it Hill-type formula because a similar 
formula for periodic orbits was shown by Hill when he considered the motion of lunar perigee IH at 1877. 
However, Hill did not prove the convergence of the infinite determinant, and the convergence was given by 
Poincare ifTSl . Afterwards, the Hill-type formula for a periodic solution of Lagrangian system on manifold 
was given by Bolotin (ll. For more results, please refer tll2ll.ll4ll. lfT0ll .ll3l') etc. We should point out that, till 
now, all the known results on the Hill-type formula were given for the periodic-type boundary problem. 

To state Hill-type formula for the separated boundary conditions, we firstly introduce some notations. 
Suppose A is a Lagrangian subspace of (R^”, a»o), a Lagrangian frame for A is a linear map Z : R.” —> 

whose image is A. It is easy to see that the frame is of the form Z = ^ j, where X,Y sue nxn matrices 

and satisfied X^Y - Y^X. 

By fhe standard Legendre transformation, the linear system ( 11 . 11 ) with the boundary conditions ( 11 . 21 ) 
corresponds to the linear Hamiltonian system. 


z ^ JB^iOz, z(0) € Ao, z{T) € Ai, 


(1.3) 


with 


BaU) = 


p-\t) 


Q{tfP-\t)Q{t)-R{t)-AR,{t) 


(1.4) 


Without confusion, for Lagrangian system, denote 7 ^(t) the fundamental solution of (11.31) . that is yxit) = 
JBA{t)yA{t) with 7i(0) = hn- Let Zo,Zi be frames of Aq, Ai. Obviously, yA{T)ZQ are frames of 7 i(r)Ao. 
( 7 ^(r)Zo,Zi) are 2n x 2n matrices. 

To simplify the notation, set A = -^,{Pjf + Q) + Q^^i+R, which is a self-adjoint operator on 
L^([0, r],R”) with domain: 


D(Ao,Ai) - {ye r],R”),z(0) e Ao,z{T) e Ai). 


Throughout of the paper, without loss of generality, we will assume A is nondegenerate, that is, 0 is not 
an eigenvalue of (IL1HL2I) . It is obvious that T is a nonzero eigenvalue of the svstem (ILlHL2l) if and only 
if is an eigenvalue of /?iA“^. In what follows, the multiplicity of an eigenvalue Aj means the algebraic 
multiplicity of /?iA“^ at -^lAj. 

Theorem 1.1, Under the nondegenerate assumption, we have 

]^(1 - Af) = det( 7 i(r)Zo,Z,) • det(7o(r)Zo,Zi)-\ (L5) 

j 

where the left infinite product takes on the eigenvalues Aj counting the multiplicity. 

Remark 1.2. The Hill-type formula for S -periodic orbits was built up in )[7j/ with the following form 

f](l - Af) - det( 7 i(r) - Sd) • det( 7 o(T) - Sd), (L6) 

j 

where Sd - diag{S,S). Although ( li.5D is similar to ( 17.6D . the proof is different. ( 17.61 ) is derived from the 
Hill-type formula for S-periodic orbits of Hamiltonian systems jj^. The proof of M.5^ is direct, and could 
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cover the case of di.6l) . To the best of our knowledge, di.5D is the first study on the Hill-type formula of 
non-periodic type boundary problem. The corresponding formula in Hamiltonian systems for the orbits 
with Lagrangian boundary conditions is still open. 

Trace formula is a powerful tool in the study of eigenvalue problem, especially in estimating the first 
eigenvalue. The first work on the trace formula was established by Krein lfT^ fldll for the -1-periodic orbits 
in the simple case. For the system with S -periodic boundary condition, the trace formula was established 
in Eiini. The present paper is a continuous work of Q [T^, and we will build up the trace formula for 
separated boundary value problem of Sturm-Liouville system. The idea to get the trace formula is similar 
to that in 171. That is, using T/?! instead of /?i, and give Taylor expansion on both sides of the Hill-type 
formula. With the notations defined in Section [3l we have the following theorem. 

Theorem 1.3. Assume A is non-degenerate, Aj are eigenvalues of the Sturm-Liouville system f li.iHi.2l) count¬ 
ing the multiplicity, we have for any positive integer m, 

= Z rKG„--Gj.)]. (1,7) 

j j k=l j[+-+ji,=m 

where Gk, k e N, defined in d3.4l) . are nxn matrices. 

For applications, a main observation is that the trace formula can be used to estimate the non-degeneracy 
of the system. Moreover, we can estimate the relative Morse index and Maslov-type index. The Maslov-type 
index is a powerful tool in study the stability problem, please refer IIT2l . lfm .ll8l for the detail. By using the 
trace formula and Maslov-type index theory, in Q we studied the stability region and hyperbolic region of 
elliptic Lagrangian orbits in planar three body problem. 

It is not hard to see that all the results on the applications of trace formula in Q have twins in the case 
of separated boundary conditions. We will not list all the theorems here, but give some computations for 
the trace formula some special case. Let 7? be a continuous path oi nxn symmetric matrices on [0, T], we 
consider the system 

y -t d/?y = 0. (1.8) 

Set /?"*■ = i(/? -I- 1/?|), which is a path of nonnegative symmetric matrices, we have the following corollary. 
Corollary 1.4. Suppose Tr(^ < 1, the Dirichlet problem for di.SD has no nontrivial solution. 

This result can be used to estimate the conjugate point alone a geodesics in Riemannian manifold. For 
reader’s convenience, we give details here. Let c : [0, a] ^ Af be a geodesic of Riemannian manifold 
Af. Choose {ei(0), ...,c„(0)} to be an orthogonal normal basis of cCO)-*- c rc(o)Af. Its parallel transport 
{e\{t), ..., en{t)}{t e [0, a]) along c gives an orthogonal normal basis of c{t)-‘- c Tc{t)M. Recall that the Jacobi 
equation is 

— - + R{c, J{t))c = 0. 
dt 

The point cff) is said to be conjugate to c(0) along c, to ^ [0, a], if there exists a nontrivial Jacobi field J 
along c, with 7(0) = /(to) = 0. Suppose /(t) = YIi=\ 7((t)c,(t) is the Jacobi field along c, fhen Jacobi equation 
can be rewritten as 

n 

fit) + J] Rij{t)Jj{t) = 0, / = 1,..., n, (1.9) 

t=i 
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where Rij{t) = {R{c{t),ej{t))c{t),ei{t)). Let R{t) = {Rij{t)), which is a symmetric matrix, then c(to) is 
conjugate point if and only if the second order system X{t) + R{t)X{t) with Dirichlet boundary conditions 
has a nontrivial solution on ? € [0, tol- It is obvious that R{t) := Tr{R{t)) is the Ricci curvature in the 
direction of c. Set R^{t) = \{R + |R|), then Corollary 1 1.41 implies that there is no conjugate point alone [0, a] 

if (L 

Now, if we consider the system (11.81) in the case n = \,R- 1 with the boundary conditions 

y(0) - 0, cos(%(r) + sin(0)y(r) = 0, 0 e [ 0 , 7 r/ 2 ]. (1.10) 

It is well known that the k-th eigenvalue Ak is the k-th positive solution of the next transcendental equation 

tan( Vdr) = - tan( 0 ) Vd. 

It is easy to check that if 0 = 0, then Ak = and \f 6 = njl, then Ak = ^{k - For 9 € (0, 7 r/ 2 ), it is 
obvious that {k - < ^fAkT < kn. However, Ak can only be solved numerically. As an application of the 

trace formula, we have the following equality, which itself is interesting. 

y 1 sin(g) + cos(g) 

^Ak 6 (sin( 0 ) + r cos( 6 »)) 

Obviously, for 9 = 0, (11.111) gives the well known identity Z 7 ^ - and for 9 = nil, (11.111) gives the 

jteN ^ 

identity ^ ^ ^ , = i. To the best of our knowledge, for 9 e (0, n'/2), we don’t know any such kind of 

keH 

formula before on the sum of The detailed calculation will be listed in Section 4. Moreover, it is worth 
to point out that we can compute the value of X for any ni € N by the trace formula (11.71) . 

''■k _ 

The present paper is organized as follows. In section 2, we give the proof of the Hill-type formula (11.51 ). 
Section 3 is devoted to proving the trace formula. Finally, in Section 4, we will give the proof of Corollary 
1 1.41 and identities (11.111) . 

2 Hill-type formula for Sturm-Liouville systems 

In this section, we will give the proof of the Hill-type formula. The following lemma coming from |[T^ 
Lemma 3.6] plays a important role. 

Lemma 2.1. Let f{z) be an entire function with zeros at z\,Z 2 , ■ ■ ■ (counting multiplicity). Suppose f satis¬ 
fied 

i) Exponential bounded condition: for any e, there exist such that 

\f{z)\ < Cfexp(elzl), 

ii) Sum finite condition: then 

00 

/(z) = /(O) ]~[(1 - z~^z). 

1=1 
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Since is a trace class operator, by ifT^ Chapter 3, P33], we have that det(/ + is an entire 

function with |det(/ + /t/?iA“')| < exp(|/l| • ||7?iA“^||i), where || • ||i is the trace norm. Moreover, it satisfied 
the exponential bounded condition. It is obvious that A,, is a zero point of det(/ + if and only if A„ 

is an eigenvalue of the system (I1.1H1.2I) . It follows that Z iin < where the sum takes for Aj counting 

n 

multiplicity. From Lemma l2Al we have that 

det(/ + d/?iA-') - ]~[(1 - (2.1) 

n 

To continue, it is easy to verify that yo is a solution of (I1.1H1.2I) with respect to the eigenvalue To if and 
only if zo is a solution of (11.31) with respect to the same eigenvalue. And it is equivalent to yAQiT)zoiO) € Ai. 
We have the following observation. 

Lemma 2.2. Suppose that A is nondegenerate, then 

dim ker(/?iA“' + I/Tq) = di}nyAf,iT)Ao n Ai. (2.2) 

For the Lagrangian frames Z, of A,, / = 0,1, set 

g(T) = det(r,(r)Zo,Zi). (2.3) 


We have 

Lemma 2.3. g{A) is an analytic function and satisfied the exponential bounded condition. 

Proof. The analyticity of y^ comes from Krein iflAll essentially. For the Taylor expansion, readers are 
referred to I?) Section 2.2]. Next, we will show that g(A) satisfies the exponential bounded condition. For 
nonzero A, let p = A^^^ and a{p) - diag{pl„,p~^ln), we set y,i(f) = a{p)~^yA{t), direct computation shows 
that 

TSy^iT)) = Ja{p)BA{t)a{p)y{t)- 

dt 

Moreover, 

a{p)BA{t)a{p) ^ Bf,{t) + p^B{t) 

with 

0„ -P-'G \ 0« 

_QTp-l p-\QTp-^Q-R) y 0„ -Pi 

Let y^ be the fundamental solution with respect to then 

= p^Jylht)y^ ■ y^jAit)). 

Restricting on the region |/i| ^ 1, it is obvious that is bounded, and thus yJB(t)y^ is bounded. Hence 
||y“^y^(r))|| ^ exp(C|/ip) for some constant C. Consequently 

llyi(T)II^Co|T|i/2exp(C|T|i/2). 

Finally, notice that g(A) is the finite combination of the finite product of the branches in the matrix, we have 
the results. 

□ 
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By Lemma \1J[ ^(Aq) = 0 if and only if Aq is eigenvalue of (I1.1H1.2I) . That is, g(A) has the same zero 
points as det(/ + /l/?iA“*). Moreover, we have the following lemma. 

Lemma 2.4. Suppose Ri > 0 and Aq is a zero point of g(A), then the multiplicity ofg{A) at Aq is same as the 
multiplicity o/det(/ + at Aq. 

Proof. Suppose the multiplicity of g(/l) and det(/ + /l/?iA“‘) at do is m\ and m 2 respectively. Since R\ > 0, 
the eigenvalue of is simple and then m 2 - dimke.x{I + By Lemma [T2l we have m 2 ^ m\. 

On the other hand, by the techniques of small perturbation (details could be found in ||9l Section 4]), we can 
assume g(A) has mi simple zeros near do, and thus mi ^ m 2 , which implies the result. □ 

From the above lemmas, we will give the proof of Theorem 1.1. 

Proof of Theorem [77/] We firstly prove the Hill-type formula for the case Ri > 0, by the nondegenerate as¬ 
sumption, 0 is not a zero point of g(A). Please note that both det(/ - 1 - AR 1 A“*) and g(d) satisfy the exponential 
bounded conditions and by Lemma [2^ they have the same zero points with same multiplicities. Next by 
Lemma [2A1 we have 


det(/ + d/?iA-')-g(0)-ig(d). (2.4) 

In the general case, choose oro e R such that Ri - a^ln > 0 and A - 1 - a^ln is nondegenerate, then 

det(/ -t d/?iA“‘) - det[/ -t A{Ri - aoIn)(.A + aoInT^] ■ det(/ -h (2.5) 

By using (12.41) on the two factors of the right hand side of (12.51) . we have (12.41) for the general /?i. By taking 
d = 1 we get the desired result (11.51) . □ 

3 Trace formula for Sturm-Liouville systems 

In this section, we will prove Theorem 11.31 To do this, we will consider the expansion of the Hill-type 
formula (12.41) . Notice that A~^ is a trace class operator, by ifT^ P47, (5.12)], 

det(/ -t d/?iA“‘) - exp( — - -d“rr((/?iA“')'")). (3.1) 

m=l 

Next, we will give the expansion on g(d). Let Vq = Aq n Ai, assume dim Pq - ^o^ then 0 < ko < n. 
Suppose that {di, ■ ■ ■ ,(7^^) is an orthonormal basis of Vq, and {di, ■ ■ ■ ,dk^^,dkQ+\, ■ ■ ■ ,dn} is an orthonormal 
basis of Aq. Notice that R^” = Aq 0 /Aq. Therefore, setting d^+j = Jdj, we have {di, - ■ ■ ,d 2 n} is a basis 
of R^” and the matrix Mi = {di,d 2 , - ■ ■ , d 2 n) is a symplectic orthogonal matrix. Next, set Pi = Ai 0 Pq, 
then it is a Lagrangian subspace of R^” 0 (Po 0 dPo)- Take an orthonormal basis {fk^+i, • • • ,/«) of Pi, then 
{di,--- ,dk^, fkf^+i, ■ ■ ■ , /„) is an orthonormal basis of Ai. 

Let [ek',k = 1, • • • ,2?i) be the standard basis of (R^”,mo). Obviously e„+k = and M^dj - ej for 
I ^ j £ n. Notice that M^{fkg+i, ■ ■ ■ ,f„) gives a Lagrangian frame of M^Vi. By direct computation, for 
Icq + I S I ^n and 1 ^ 7 ^ ko, 

{M\fuej) = {Mlfue„^j) = 0. 
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Rewrite such a frame as 
/-I 


M 2 = 


^n-ko ^n-ko 


Xi 

Yi 

, and 


, where Xi, are (n - ko) x {n - ko) matrices and Fi is nonsingular. Let 


Xfs - ihko XI 2 ) ■ Ml, 


(3.2) 


^ 4 „ 0 0 0 ^ 


where l 2 ko o M 2 = 


. Obviously, M 3 is a symplectic orthogonal matrix. 


0 In-ko 0 -liF -1 

0 0 4q 0 

0 0 0 In-ko 

Let Vo = span{ei,- ■ ■ ,c„}, Vi = span{eko+i, ■ ■ ■ ,eko+n} which are Lagrangian subspaces. Let Pq, Pi be 
the orthogonal projections onto Vo and Vi respectively. For any matrix M on we always set 

P(M) := P 1 M 3 MM 3 ‘Po, 

which is a n X n matrix with P{M)i j = {M^MM^^ej,ei+ko). In the case dim Vo = 0 or n, the expression 
of P(M) is simple. In fact, rewrite M^MM^^ - 


/A A •. 

^ Ml M 2 ^ 


M 3 M 4 


, then P(M) = M 3 in the case Ao = Ai, and 


P(M) = Ml in the transversal case Ao n Ai = (0). 
Notice that det(M 3 ) = 1, 


g(d) = det(M3)det(r2(P)Zo,Zi) - det(M3ri(P)Zo, M3Z1). 

Direct computation shows that 

det(M3r2(r)Zo,M3Zi) - det(M3ri(P)M3-iM3Zo,M3Zi) - (- 1 )”^« det(P(y2(P)))det(Fi). 

Then 

g(d)g(0)-i = detiPiyAT)) • det(P(yo(P))“^). 

Set D = diag{0n, -Pi), then = Po + from Q, let 


(3.3) 


D{t) = yl{t)D{t)yo{t), 


and 


’4-1 


pi ^ p‘l ^ pH 

Fk = I JD{t\) I JD{t2) ■ ■ ■ I JD{tk)dtk ■ ■ ■ dt2dt\,k € ^ 
Jo Jo Jo 

By Taylor’s formula. 


ym - yo(P)(4« +AFi + --- + A^Fk + ---), 


then 


njAiT)) = PiyoiT)) + AP{yo{T)Fi) + ■■■ + A’^P{yo{T)Fk) + •••), 
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where PiyoiT)) is nonsingular. Set 


Gk = P{yo{T)Fk) ■ nyo{T))-\ for ^ € N, 


(3.4) 


and let f{A) - det(/„ + AG\ + ■■■), which is an analytic function of A. It is obvious that g(/l)g(0)“^ = /(A). 

Since f{A) vanishes nowhere near 0, we can write f{A) = then by l|71 Formula (2.6)] and some 
direct computation, 


Jm), 


( 0 )/„,! . Tr{Gj, ■ ■ ■ Gj,)). 


(3.5) 


Compare the coefficients in (13.11 ) with (13.51 ). we have 


Tr{{R,A-Y) = Tr{Gj, ■ --Gj,)). 


k+m 


(3.6) 


k=\ 


Jl+-+Jk=m 


{-!)'" 


This proves Theorem [O] because Tr{{RiA ^)"‘) = Yj 

j ^ 

Moreover, for the first two terms, we can write it more precisely. 

-T 


Zj: = -^^(<^i) - -Tr(F(yo(T) ■ J y^(t)D(t)yo(t)dt) ■ P(yo(T))-'), 


(3.7) 


and 




- Tr{G\)-lTr{G 2 ) 

('P(yo(T)-jJ y^(t)D(t)yo(t)jJ y^(s)D(s)yo(s)dsdt) ■ nyoiT))-^} 

pT 2 

([^(ro(7’) • J y^(t)D(r)yo(Odt) • nyoiT))-^] )■ 


- - 2 rr( 
+Tr( 


(3.8) 


4 Examples 

In this section, we will give some detailed calculation on the trace formula for some special separated 
boundary value problems for Sturm-Liouville system. At first, we will consider the Dirichlet problem for 


the system (11.81) . Obviously, in this case A = R\ - -R- Let K„ = | *“ I, D = 

0„ 0„ I [ 0„ R 


7m Om 


o„ 0 „ 


. Recall 


that yo(0 satisfied yo(t) - JK„yo{t) wifh yo(0) - hn- Direcf compufafion shows fhaf yo(t) = 
is easy fo verify 7o(0~^ = I ?" I- Therefore, 


Ifi On 

tlyi If, 


.If 


-tl,i If, 


JD{t) = yl\t)JD{t)ya{t) = 


-tR{t) -R{t) 
fR{t) tR{t) 


-L^Rdt] 

Ddt= 

Jo V C fRdt f tRdt J 


Then 






and 




^‘"Kdl-TgtRd, gtRdl-rgRdt 


yo{T)-J I Ddt = 

Obviously, PiyoiT)) = Tin, and 

r(yQ{T)J bdi) = bRdt - T ^ iRdl, 


thus 


We have 


1 

Gi = — I bRdt - I tRdt. 

T Jo Jo 

Tr(RA-‘)='^± = Tr(£(,-^)Rd,). 

Recall that R^ = j(R + |/?|) is nonnegative matrices, let /f!" be the y-th eigenvalue of y + AR'^'y = 0 under 
the Dirichlet boundary conditions, then Aj > 0 for j € N. Similar to the discussion of E Theorem 4.12], 
Y,j j- = dt^ < 1 implies Ti > 1, hence ^ is nondegenerate for A € [0,1]. Since 

R^ > R, we have ^ -\-Ri^ nondegenerate. This proves Corollary 11.41 


At the end of this paper, we will consider (11.81 ) with the boundary condition (11.101) . We choose d\ = 


and /2 ^ 


cos( 0 ) 
- sin( 0 ) 


to be the frame of Aq and Ai separately. Then M\ = I 2 , M 2 = 


and consequently M 3 = M 2 . It is not hard to see, M^ ^ = 
Direct computation shows that in this case 


1 - cot( 0 ) 

0 1 


Rewrite M = 


1 
0 
a b 


!P(M) = a + ccot(0). 

So we have !P(yo(T)) = 1 + T cot(0), and easy computations show that 


ni Y'Z 'J'* 

p{n(TV I Dd,) = -- - - 


COt(0). 


We get 


G, - - 


3T^ + com 372 sin(6») + cos(0) 


6(1 + Tcotm 


6(sin(0) + T cos(0)) 


cot(0) 

1 


in short. 


(4.1) 


(4.2) 


(4.3) 


By (13.71) . we get (11.111) . It should be pointed out that maybe the identity (11.111) could obtained by some other 
method. However, by the trace formula, we can get many other interesting identities directly if we consider 
different boundary conditions for the Sturm-Liouville system. 
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